Gap theory of rectification in ballistic three-terminal conductors 
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We introduce a model for rectification in three-terminal ballistic conductors, where the central 
connecting node is modeled as a chaotic cavity. For bias voltages comparable to the Fermi energy, 
a strong nonlinearity is created by the opening of a gap in the transport window. Both noninter- 
acting cavity electrons at arbitrary temperature as well as the hot electron regime are considered. 
Charging effects are treated within the transmission formalism using a self-consistent analysis. The 
conductance of the third lead in a voltage probe configuration is varied to also model inelastic ef- 
fects. We find that the basic transport features are insensitive to all of these changes, indicating 
that the nonlinearity is robust and well suited to applications such as current rectification in ballistic 
systems. Our findings are in broad agreement with several recent experiments. 

PACS numbers: 73.40.Ei, 73.23.Ad, 73.63.Rt 
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There has been recent interest in applying ballistic 
three-terminal junctions as voltage rectifiers or diodes 
in emerging nanoelectric technology. From an appli- 
cations point of view, the observed nonlinear I — V 
(nu've-——— —————— in these junctions are quite at- 
tractive in that the effect persists to room temperature 
and originates without any special engineering. Initial 
experimental investigation into the issue sought to add 
the source of nonlinearity in various ways: The experi- 
ment Ref. [Hi placed a triangular obstacle in a cavity in 
an effort to "force" the electrons in one direction or an- 
other. For some time it was thought that such enhance- 
ments were necessarj&i^ but other discussions^ pointed 
to a more general origin of the rectification effects* 4 These 
discussions ^ 13 ' 14 are all based on an application of the 
transmission approach for multi-terminal conductors^ 5 - in 
the non-linear regime. However, beyond the linear regime 
application of the transmission approach requires a self- 
consistent treatment! 16 ' 17 Indeed, it has been shown that 
ballistic junctions are very sensitive to side-gating. 18 

Even in the linear regime ballistic four-probe junctions 
have found applications in Hall micromagentometr y 19 ' 20 
and scanning Hall probe microscopy. 21 ' 22 Recent ad- 
vances include a vector Hall sensor^ Multi-terminal 
ballistic junctions are also found to be sensitive 
potentiometers^ 

In mesoscopic physics, the properties of ballistic four- 
probe junctions were originally investigated in the lin- 
ear transport regime starting with work by Roukes et 
ali 2 ^ who found at low temperatures an absence (quench- 
ing) of the Hall effect. Different geometries were inves- 
tigated leading to an enhanced or suppressed Hall effect 
depending only on the geometry of the Hall cross. 26 ' 27 
Already at He temperatures these effects can be well 
described by classical trajectories.— Interference effects 
play a role at much lower temperatures^ 9 - In the non- 
linear regime, interference effects in chaotic cavities have 
recently found interest in connection with the generation 
of rectification effects. These works examine (predom- 
inantly) the second order in voltage term of the I-V- 



characteristioii and demonstrate that interactions lead 
to deviations from the Onsager symmetry. 30 ' 31 ' 32 ' 33 ' 34 ' 35 
In comparison, the somewhat extreme conditions of large 
bias and high temperature envisioned for applications of 
ballistic structures as rectifiers and diodes are outside 
the scope of the mesoscopic physics literature and thus 
require a separate treatment. 

The purpose of this article is to present a simple model 
of classical rectification in ballistic chaotic cavities. We 
here take a minimalist approach to the problem, and 
make only the following two assumptions: (i) Trans- 
port between lead and cavity via the quantum point 
contacts (QPC) is ballistic so the Landauer formula ap- 
plies, and (ii) The mean level spacing in the cavity is 
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1: Energy landscape for the low (A) and high (B) bias 
case. In case (A) transport is approximately linear, and the 
current-carrying electrons have energies in the transport win- 
dow. In case (B) the electrical bias V is sufficiently large that 
an energy gap opens between the transport window and the 
filled Fermi sea of the right lead. The energy gap is responsible 
for the strong nonlinearity in the transport characteristics. 
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much smaller than the charging energy of the cavity, 
and therefore charge neutrality of the cavity under non- 
equilibrium conditions is imposed. From these two as- 
sumptions we develop our model, and demonstrate that 
when the applied voltage is comparable to the Fermi en- 
ergy, a strong nonlinearity develops. Furthermore, we 
demonstrate that this behavior is insensitive to the de- 
tails of the model, showing that this mechanism is generic 
and robust. 



I. MINIMAL MODEL 

We now describe the model in detail. A ballistic cavity 
is connected via three QPCs to bulk leads. We assume 
that the energy barriers of the QPCs (see Fig. 1A) are 
specified by a potential energy Vq. Electrical bias V is 
applied across two of the leads. When the third lead is in 
the voltage probe configuration, there are two basic DC- 
transport characteristics: the dependence of the potential 
V p on the third lead (voltage probe), as a function of V, 
and the I—V curve through the left and right leads. When 
the third lead is electrically fixed to be the mid-point 
voltage between the left and right lead, there are then 
the three currents through the leads, linked by current 
conservation. 

In the voltage probe situation, an even simpler model 
that captures the basic physics is to truncate the third 
lead by pinching off the third contact, and to consider 
the dependence of the internal cavity potential U versus 
V. We will first work out this simplest case analytically 
in detail, and present only the numerical results for more 
realistic extensions of this model. 

For a chaotic cavity the occupation function is 
isotropi o 36 ' 37 and the main theoretical task is to find its 
dependence on energy. The use of a chaotic cavity gives 
the results of our minimal model a degree of universality 
which is absent in ballistic junctions with short geometry- 
dependent trajectories. 

Low bias limit, V < Ep — Vq. — We first assume zero 
temperature and elastic scattering. The current going 
into the cavity from both leads is given by the Landauer 
formula 

I a = J dEj a (E) = (e/h) J dET a (E)(f a - f c ), (1) 

where a = L,R, j(E) is the energy-resolved current, 
T(E) is the energy-dependent transmission, and Jl,r,c 
are the occupation functions of each region. We assume 
that the lead occupation functions are completely speci- 
fied by Fermi functions with a single potential: Ep + V 
and Ep. The large energy scales involved imply that the 
transmission of the QPCs may be treated semi-classically, 
so only the coarse energy-dependence is kept: 

T a (E)= (E ~ Ea) Q(E-E a ), (2) 
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FIG. 2: (A) Current / between left and right lead is plot- 
ted versus applied bias V in units of Ef for different tem- 
peratures. We take Vo = -2Ep, so the cross-over is around 
Ef — Vq = .8Ef- Increasing temperature raises the inter- 
nal potential U for large applied bias. (B) Cavity potential 
U — V/2 plotted versus applied bias V in units of Ef, for dif- 
ferent temperatures. Inset of (B): Considered geometry. Ex- 
ternal bias is applied across left and right lead, /il = fiR + eV. 
The cavity voltage is measured with the help of a third voltage 
probe. 

where we use a linear interpolation (valid for energies 
larger than the conductance quantization scale) and 
El = V + Vq, Er = U + Vq are the minimal ener- 
gies required for carriers to pass through the QPCs (see 
also Fig. IB). The energy scales A Q characterize how 
open the contacts are. The limit A — > oo corresponds 
to a closed contact. Other energy dependence (such as 
the semi-classical 3/2 transmission law) leads to similar 
physics. 

Imposing energy-resolved current conservation ji,(E) + 
Jr(E) = 0, allows us to solve for the cavity occupation 
as a function of energy: 

E-V-Vp 
/low ~ 2E-2V -V-U' ( ) 

where Ep < E < Ep + V. Between U and Ep, the cavity 
occupation is 1. We can now calculate the charge in the 
cavity, 




where D is the (constant) density of states, V a is the 
potential of lead a, and C' a is the capacitance linking the 
central cavity to terminal a. For simplicity, we focus on 
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the realistic case of reasonably large cavities, where the 
cavity mean level spacing, D , is much smaller than the 
charging energy of the cavity, e 2 /C. In this case, Eq. (j4|) 
is a charge neutrality condition, and the left hand side 
may be replaced by the equilibrium charge on the cavity, 
Qo = eDEp. Inserting Eq. © into (gj), we find 

Q/(eD) = E F = E F - U + V/2 (5) 

(2E F -V -U -2V 



l(U-V)/4] log 



\2E F - U + V-2Vo 



giving a self-consistent (transcendental) equation for the 
unknown potential U, that may be solved numerically. 

Notice that the logarithm in Eq. ([5|) originates from 
keeping the energy dependence of the transmission, and 
integrating over the full energy range. Indeed, if we sim- 
ply neglect the logarithmic correction, we recover the 
usual linear result, U — V/2. In Fig. 2(B) we plot this 
solution for V < Ep — Vq (the T = curve only). In 
order to take out the usual linear behavior and focus on 
the nonlinearity, we plot U — V/2 (all other plots of cavity 
or probe potential will reflect this convention). 

High bias limit, V > Ef — Vq. — In this range of param- 
eters, the large applied bias lifts the energy barrier of the 
left QPC up above the Fermi level of the right contact 
(see Fig. 1(B)). This opens up a gap in the transport 
window, corresponding to an energy range too low for 
the left current carriers to fill. Here 



/high 



U <E < E F , 

E F < E < V + V , 

V + V < E < E F + V. 



(6) 



Inserting this to Eq. ([4]), we find the equation 

Ep = Ep — U - 
[(U-V)/4] log 



Q/(eD) = Ep=E F -U+(Ep- V )/2 (7) 

2E F - V -U -2Vo 



2E F -U + V-2V 



We observe that if the logarithmic contribution in Eq. 
(7) is neglected, we obtain U = (E F — V )/2, so U goes to 
a constant as V continues to increase. This explains why 
Fig. 3(A), (C) as well as the experimental data 11 (which 
plot U — V/2) show transitions from a flat dependence to 
a shifted line with slope -1/2 at large bias V . The energy 
scale E F — Vq is the cross-over point. The numerical 
solution of U as a function of V is given in Fig. 2(B) for 
V > E F — Vq (T = curve only). Note also the spatial 
inversion symmetry, V — > —V,x — > —x, for a symmetric 
geometry gives the (trivially) symmetric negative voltage 
behavior as found in recent experiments^ For GaAs, the 
Fermi energy is around Ep ~ 20 meV, which is of the 
same order of magnitude as room temperature. 

Finite temperature. — At finite temperature, again im- 
posing energy-resolved current conservation and solving 
for the cavity occupation, we find the well-known result 



but this only applies for E > V + Vq ■ For lower energies, 
current can flow only from the right lead (or not at all), 
so here fc = fn, the equilibrium Fermi function. 

Calculating the total charge in the cavity as before, we 
find 



rVo + V pc 

Q/(eD) = dEf R + 

JU JVn 



IlTl + JrTr 



-v +v 

dE- 

tu JVo+v Tl + Tp;, 

Tlog[l + exp(E F /T)}. (9) 



This gives a self-consistent equation for the cavity poten- 
tial U, as a function of V, T. The results are shown in 
Fig. 2(B) for different values of temperature T. Higher 
temperature tends to increase the cavity potential. 

It is also interesting to look at the I— V curve for trans- 
port between the left and right lead. From Eq. (fTJ) , we 
can now find the current as a function of applied voltage, 
now that we know the internal potential U: 



e f°° 

- dET L (f L -f c 

n JV+Va 



'V+Vo 

h 



dE 

V+Vo T L + Tr 



Ul - fn), 



(10) 
(11) 



where the transmission Tr of Eq. ^ depends of the cav- 
ity potential U. Fig. 2(A) shows the I — V curve as 
temperature is varied. Higher temperature tends to in- 
crease current and smooth the transition. The current is 
an antisymmetric function under voltage reversal. 
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FIG. 3: Rectification characteristics as the width of the volt- 
age probe (A/Ap) is varied for T = (A), (B), and T = .5E F 

(C) , (D). As in Fig. 2, Vb = -2Ef, so the cross-over is around 
Ef — Vo = .8Ef- As the temperature increases, the cross-over 
is more smooth for both the probe voltage (C) and the current 

(D) . When the probe is turned on from closed (A/A p = 0) to 
open (A/Ap = 1), the T — curves are essentially indifferent, 
while the T — .5Ef curves have slightly more variation. 



fc = (T L h + T R f R )/(T L + T R ), 



(8) 
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II. VOLTAGE PROBE MODEL 

We now open up the third lead on the structure. 
Adding the third lead does two things: First, it allows a 
realistic probe of the cavity voltage, which is the essence 
of the rectification effect, and second it allows a way of 
treating inelastic processes in a phenomenological way. 
This latter effect occurs because high energy carriers can 
enter the probe, and be reinjected at lower energy^ In 
reality, charge redistribution in energy occurs not only 
due to the voltage probe but through various inelastic 
scattering mechanisms i 39 ' 40 In our model the probe can 
be turned on and off with a coupling parameter. 

The equations are similar to those of the previous 
section: Now there are three energy-resolved currents, 
3LijR)3p where j p is the energy resolved current enter- 
ing the voltage probe. We require 

• (i) Energy-resolved current conservation, +jn + 

Jp = o. 

• (ii) Charge neutrality, Q — eD J dE fc(E) — 
eDE F . 

• (Hi) A new condition: No net current flow in or out 
of the voltage probe, J dE j p (E) = 0. 

These equations determine both cavity and probe volt- 
age. We take j p (E) = (e/h)[(E -Vq- V p )/A p ](f p - f c ). 
The conductance of the probe lead can be smoothly 
turned on and off by varying A/A p from (the off con- 
figuration), up to oo where it dominates over the other 
leads. 

Low bias limit, V < Ep — Vq. — Solving for the cavity 
occupation from (i), we find 



fc 



1 

fc 



U < E < Ep, 
Ep <c E < Ep 



v, 



v 



(12) 



where 
fc = 
ft = 



f+ E F + V p <E<Ep + V, 
E — Vq — V + (A/A P )(E -Va- V p 



2E-2V -V-U 
E- 



h(A/A p )(E-V 
V Q -V 



v P ) 



2E - 2V - V - U + (A/A P )(E - Vq 
The first equation from (ii) is given by 



V P ] 



£13) 



rE F + V p pHip+V 

Q/eD = E F -U+ dEf c + dE f~Q = Ep. 

JE F JE f +V p 

The second equation from (Hi) is then 



E F + V t 

E F + V P 



V +V p 



dE (E — V — V P )(E -Vp-V) 
2E - 2V - V - U + {A/A P )(E -V - V p ) 



dE(E-V -V p )(E-V -U) 



2E -2V -V 



U+(A/A P )(E-Va-V p )' 
(15) 



As before, these integrals may be expressed as loga- 
rithms, but the equations are again transcendental, so 
they need to be solved numerically. 

High bias limit, V > Ep — Vq. — Here also, the main 
feature is the presence of a gap in the transport win- 
dow. There are four energy windows, two of which are 
nontrivial: E € [U, Ep, Ep + V p , V + V , Ep + V]. 

• Filled region: Between E G [U, Ep] every state is 
occupied, so fc = 1. 

• Low energy region: In the range E € [Ep, Ep+V p ], 
the voltage probe can inject carriers into the cavity. 
Solving for the occupation in this range from energy 
resolved current conservation, we find, 



f LE 

JC 



E-V p -V 



E-V p -Vq + (A P /A)(E -U-V ) 



(16) 



• Gap region: Between E G [Ep + V p , V + Vq], the 
left electrons are injected at too high an energy 
to fill this region, so jh = 0, while the energy of 
the right and probe electron is too low to fill it: 
JR = (e/h)[(E -U- V o )/A](0 - f c ), and j p = 
(e ■■/ h)[(E - V p - V ) I 'A p ](0- f c ). Therefore f c - 0, 
giving the gap. 

• High energy region: In the range E G [V + Vq, V + 
Ep], there is injection from the left, and drain to 
the right. Current conservation yields 

f HE = E-V-Vq 

Jc 2E-2V -V-U+(A/Ap)(E-V p -V Q )' 

(17) 

Notice A, A p have switched places from the low 
energy region. 

The first equation from (ii) is given by 

Q/eD = E F = E F -U (18) 



V P +E F 



dEf 



LE 







E F +V 



dEf 



HE 

C ■ 



lEw JVo+V 

The second equation from (Hi) is given by 

[ Ef+v * dE (E - V p - V a )(E - Vp - U) 
Je f E-Vo-U+(A/A p )(E-V -V p ) 



(19) 



Ep+Vr, 



dE (E-V p -V )(E-Vo-V) 



Vo+Vp 2E - 2V a - V - U + (A/A P )(E -V a -V p \ 

(20) 

Finite temperature. — The analysis is somewhat simpler 
at finite temperature, simply because there are fewer en- 
ergy regions to keep track of. Now there are three energy 
regions to attend to: E G [U, Vq + V p , Vq + V, oo] . Current 
conservation gives three possible answers for the cavity 
occupation fc, depending on what energy region we are 
considering, 




Vn 



U < E < Vq 
Vq 

Vq + V < E < oo 



V P <E<V 



V 



(21) 
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where 



/hit 
/high 



Tpfp + T R f R 
T p + T R ' 

2p/jj + Tl/l + T R j R 



T P + T L + T R 
This gives the charge in the cavity as 



(22) 
(23) 



Q/eD 



V p +V rVo+V 

dEf R + / dEf int , (24) 

Jv +V p 

/ rf^/high = Tlog[l + exp(E F /T)}. 

JVn+V 



IVo + V 

The no-net-current probe condition (Hi) then reads 

•V +v 



/fVo+V 
dET p (f p - f c ) = + / d£ T p [/ P - / int ](25) 
JVo + V p 

/>00 

+ / dET p [f p - f high ] =0. 

Taking / p to be a Fermi function with unknown energy 
V p , these two equations may be numerically solved to give 
U and V p as a function of the parameters T, V, A, A p . 

The results are given in Fig. 3. We see that the ba- 
sic rectification features remain, and changing the model 
parameters does not alter the basic picture, indicating a 
robust voltage rectification effect. 

III. CURRENT RECTIFICATION 

Current rectification occurs whenever there is a net DC 
current produced by an external AC voltage source. We 
now consider this situation in the three-terminal geome- 
try, where the AC voltage signal has a frequency slower 
than the RC time of the cavity. This situation may be 
analyzed by investigating DC transport with chemical 
potentials fiL = V/2,fi R = — V/2, /z p = 0, and how the 
probe current depends on V. When V — > — V, the left 
and right leads will switch roles, but because of the reflec- 
tion symmetry in the problem, the finite probe current 
produced by the strong nonlinearity will remain unal- 
tered. Under repeated sign changes the system will sus- 
tain a net DC-current from the probe out into the left 
and right leads. To be consistent with the previous re- 
sults, we add V/2 to all potentials, and define all currents 
as positive when they enter the cavity. 

It is instructive to estimate what is the theoretical up- 
per limit of the speed of a ballistic rectifier. This is con- 
trolled by the i?C-time (trc) of the cavity which con- 
trolles the relaxation of cavity charge: if the external 
AC frequency is slower than t R c then current rectifi- 
cation will occur because the cavity has time to estab- 
lish a (nonequilibrium) steady state, while in the oppo- 
site limit, rectification is expected to not occur. Recent 
measurements for the capacitance of a mesoscopic cav- 
ity found C ~ Taking on the order of 10 open 



channels gives a resistance R ~ lkCl. We then estimate 
the charge relaxation time to be t R c ~ 10~ 12 s. These 
parameters give rectification on a Terahertz scale. 

Making use of previous results, everything is the es- 
sentially the same, except that the probe voltage is now 
fixed, V p = V/2, and the current through the probe needs 
to be calculated. The cavity occupation fc(E) written 
in Eq. (T2"T]) is the same, as is the charge in the cavity 
Eq. (EH), but now with V p = V/2 in both (|21|24|) . These 
conditions set the potential U. We find that the current 
through the left lead is 



h 



dE 



V+V 



TL[T p (f L -f R )+T R (f L -ff. 

T L + T R + T p 



the current through the right lead is 



(26) 



/,. = £ / " dE T B{Tp{f R - fp) +T L (f R - f L )} 



V + Vo 

v+v 



dE 



T L + T R + T p 
T R T p (f R - fp) 



'Vo+V/2 T R + T p 

and the current through the probe is 



(27) 



In = 



dE 



V+Vo 
V+Vo 



dE 



T p [T L (f p -f L ) + T L (f p -f R )] 
T L + T R + T p 

T R T p (f p - f R ) 



Vo + V/2 



Tr + Tp 



(28) 



We can now use the charge equation to find the poten- 
tial U, put this into Eqs. (|26I27I28[) , and find the currents 
in all leads. The results are shown in Figs. 4 and 5 vary- 
ing temperature and the width of the probe contact. At 
small bias, the I — V curves are linear in both the left 
and right lead, with no current passing through the probe 
lead. This is the expected situation for linear transport. 
For larger voltages, as the probe lead is gradually opened, 
the current through the left lead is essentially unchanged, 
but the current through the probe and right lead both in- 
crease. This is the direct analogue of the probe voltage 
Vp and cavity voltage U following the lower right voltage 
in the probe configuration. 

Hot electron regime. — One weakness of the above anal- 
ysis is that the gap produced in the energy window is un- 
realistic at high temperatures. As electrons collide with 
one another and phonons, they will redistribute them- 
selves in energy when any inelastic processes are intro- 
duced. One essential point that must be demonstrated 
is that the nonlinearity we have discovered is not frag- 
ile to a reshuffling of electrons in the energy space. To 
this end, we will now consider the "hot-electron regime" , 
an effective model of electron transport when there is 
conservation of both charge and energy currentsj 42 ' 43 In 
the limit where electron-electron interactions are very 
strong, r ee < TCt the cavity comes to a local nonequi- 
librium steady state, described by a Fermi function with 
two parameters, [i c and T c . These parameters are deter- 
mined self-consistently by imposing current conservation 
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FIG. 4: Current rectification effects in each of the three leads for T = 0. Different curves vary the opening of the probe lead. 
For small applied bias, the I — V curves (A), (B) start out linear, while the probe lead (C) carries no current. At large bias, 
the probe lead begins to carry current (so long as it is open, A/A p > 0), and this current comes primarily by increasing the 
current carried by the right lead (the one with low-energy carriers). Current from the three leads sum to zero. Clearly, the 
influence of the probe lead plays a much stronger role here than in Fig. 3. Inset of (A): Geometry of the current-rectification 
set-up. The left, probe and right lead have applied bias [il = V, [i p = V/2, = 0. 




FIG. 5: Same as Fig. 4, for the case where T — .5Ef. Apart from smoothing the transitions a little, finite temperature also 
increases the overall scale of the effect. 
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Applied Voltage (E p /e) Applied Voltage (E p /e) Applied Voltage (E p /e) 



FIG. 6: Current rectification for the left (A), right (B), and probe (C) lead, using the hot electron model, for T = .5Ef- 
Despite the fact that the energy distribution of the cavity electrons is completely altered from the noninteracting case, the 
curves are still very similar to Fig. 5. 
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and energy conservation for transport through the cavity. 
The two constraints may be written as 

fdE^ 3a (E)=0, fdE^ Ej a {E) = 0. (29) 

Oi Oi 

Once the parameters fi c and T c are found as a func- 
tion of the tunable parameters of the system, the current 
through all leads may be found (in the current rectifier 
mode), or the voltage of the probe may be found (in the 
voltage probe mode). 

The current rectification results are shown in Fig. 6 
for the hot electron regime. We note that the trends in 
the data are the same, and the only discernible difference 
is a slight change of shape in some of the curves. This 
indicates that the model described here is robust, and not 
sensitive to changes in the cavity occupation function. 

IV. CONCLUSIONS 

We have proposed a model of ballistic rectification for 
a three-terminal geometry. The model is minimal in the 
sense that we have only included the most important ef- 
fects, and therefore should be considered a benchmark 
theory, rather than designed to predict detailed exper- 



imental features. The most important feature of our 
model is a cross-over from a weak to a strong nonlinearity 
when the bias voltage is comparable to the Fermi energy. 
This effect has already been observed in experiments. 11 
The origin of the strong nonlinearity is the opening of 
a gap in the transport window. This happens when the 
applied bias elevates the left QPC energy barrier above 
the Fermi energy of the right contact. We have demon- 
strated that varying parameters in the model, as well as 
considering inelastic effects do not alter the basic fea- 
tures of the model. Taken together, this theory indicates 
that three-terminal ballistic cavities provide robust recti- 
fication that may be used in the development of ballistic 
nonlinear elements such as rectifiers and diodes. 
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